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Abstract 

The discrete R-parity (Rp) usually imposed on the Supersymmetric (SUSY) 
models is expected to be broken at least gravitationally If the neutralino is 
a dark matter particle its decay channels into positrons, antiprotons and neu- 
trinos are severely constrained from astrophysical observations. These con- 
straints are shown to be violated even for Planck-mass-suppressed dimension- 
five interactions arising from gravitational effects. We perform a general anal- 
ysis of gravitationally induced Rp violation and identify two plausible and 
astrophysically consistent scenarios for achieving the required suppression. 



A discrete symmetry, called R-parity (Rp) is imposed on the Minimal Supersym- 
metric Standard Model (MSSM) [[I]] to ensure proton stability. This assumption 
makes the Lightest Supersymmetric Particle (LSP) stable. The most natural can- 
didate for the LSP in SUSY models is a neutralino. Indeed, calculations in SUSY 
models with soft breaking terms and radiatively induced electroweak breaking lead 
to neutralino as LSP for a wide range of allowed parameters. Moreover, the relic 
neutralino density satisfies all requirements for being the cold dark matter over 
large parameter regions. The hypothesis of neutralino as Dark Matter (DM) parti- 
cle is amenable to experimental verification []2p. Such neutralinos can be detected 
both directly through their elastic scattering off-nuclei [||] and indirectly through 
the products of neutralino annihilation [|4]],[|5p,[|6p. 

Although R-parity may remain unbroken if it is a remnant of a gauge symmetry 
[0, ||, [Tll|, [Op, there is no deep theoretical reason requiring it to be a symmetry 
of nature. In fact, many models of R P violation have been proposed [fl^l. These 
typically lead to large Rp violation and thus do not allow the neutralino to be a 
dark matter particle. In contrast, if Rp violation is very mild then the lightest neu- 
tralino is unstable, but very long-lived. Naively, one would think that a neutralino 
with a lifetime of the order of the age of the universe could be a viable dark matter 
candidate. However, most neutralino decays into visible channels, e.g. containing 
positrons, antiprotons, gamma's and neutrinos are severely constrained from ob- 
servations. They typically require neutralino lifetimes much larger than the 
age of the universe. This can only be realized if the violation of Rp is extremely 
weak []13|, [14]]. Such a scenario was already considered in the context of a specific 
mechanism [fT5|] for spontaneous R P violation driven by a vacuum expectation value 
(VEV) for the right-handed sneutrino. This mechanism requires very small Yukawa 
coupling, i.e. fine-tuning. A more natural way to obtain very small R P violation is 
to ascribe it to gravitational interactions. This violation will be described here by 
non-renormalizable terms suppressed by inverse powers of the Planck mass and 
should thus be naturally very small. 

A somewhat surprising conclusion that emerges out of our study is that dimension- 
five interactions are in conflict with the astrophysical constraints by few orders of 
magnitude if the neutralino provides the cold DM. 

Let us first quantitatively discuss constraints on the strength of R P violation. 
The lightest neutralino x is given by a superposition of wino W, bino B and two 



Higgsinos, H 1 and H 2 as: 



X = Z xWa W 3 + Z X§ B + Z^H, + Z xS2 H 2 . (1) 

We can parameterize the effective R P violating interactions responsible for the 
neutralino decay in terms of the MSSM fields as follows: 

W eff = X 1 (U C D C D C ) F + X 2 (LLE C ) F + X 3 (QD C L) F (2) 
+ e(LH 2 ) F 

The notation for the fields is standard. We suppressed the generation indices. Eq. 
@ includes only renormalizable terms relevant for neutralino decay. However, 
as it will be understood below, this expression has wider generality. 

The above interactions result in neutralino decay to three fermions. The width 
for this decay depends upon whether it proceeds through the Higgsino or gaugino 
component. In the former case, 

-i_ A 2 7 2 G F m 2 f m\ 

where m x , rhf and mj are neutralino, sfermion and fermion masses, respectively. In 
the case of quarks the width should be multiplied by the number of colours. When 
the decay proceeds through the bino component B of the neutralino the decay width 

is 2 5 

T - 1 = \ 2 7 2 aemY fR m x (A) 
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where Yj R is hyper-charge of the right fermion and 6 W is the weak mixing angle. 
Finally, the width of x — ► v + e + + e" due to the last term of Eq. @ is 

^ 1 = ^\ + ™ 2 W + I sin* e w )^£, (5) 

The condition r x > to (where to is the age of the universe) leads to constraints on Aj 
and e. However, much more stringent limits follow from production of positrons in 
our Galaxy [ |I3| , [HQ , from the diffuse gamma-radiation [ |I3| , [Hp and from neutrino- 



induced muons QT5|] . From the analysis given in QT3] ] it follows that when the decay 



to positrons is unsuppressed as in the present case, the strongest constraints on 
both Aj and e follow from the observed flux of positrons in our Galaxy. The lower 



limit on neutralino lifetime from this flux is QT4p 

t x ( x -> e+ + anything) > 7 x 10 w ^^t h, (6) 
where mi o = m x /100 GeV and h is dimensionless Hubble constant. 



Using this limit and keeping in mind indirect production of positrons through 
decay of other particles, we find 

A, < 4 x lO-^K^fti^V (— (7) 

x " \ Lev m x rrif 

e < 6 x 10- 23 Z-± m -^ 4 GeV (8) 

The above estimates show that if the neutralino is a DMP the R-parity violating 
parameters Aj and e are very strongly limited from above. 

Let us now turn to a systematic dimensional analysis of R-parity- violating op- 
erators in the MSSM. We have already discussed the operators of dimension d=4. 
Extremely small couplings are needed for the neutralino to be a DMP, which re- 
quires fine tuning. 

The relevant Planck-mass suppressed operators of dimension d=5 in the MSSM 
are the following ones: 

^{H^E^n $-(QL*U c ) D $-(U c D c *E<) D 
^(QQQH 1 ) F ^{LH^H^p . 

The first three terms in the above equation contain the F-terms of the anti-chiral 
fields H^L* and D c * respectively. These are determined in the supersymmetric 
limit by the standard MSSM R P conserving superpotential. Using this one can 
show that their contribution to the neutralino decay is suppressed either due to 
small Yukawa couplings or kinematically when the charged Higgs is heavier than 
the neutralino. The contribution of the fourth term with three Q fields to neutralino 
decay can also be shown to be small. 

The first four operators thus satisfy astrophysical constraints without need for 
extra suppression in the corresponding coefficient. In contrast, the decay induced 
by the last term cannot be suppressed kinematically. It leads to the effective R P 
breaking parameter 

e ~ f3 5 M 2 EW /M P ~ (3 5 1(T 15 GeV 

This value of e is extremely small and leads to a neutralino lifetime longer than 
the age of the universe. But it is in conflict with the astrophysical constraints @ 
unless Ps < 10~ 5 — 10 7 . This situation is similar to the gravitationally induced 
axion mass pT7|], [ |18| ] where the quantum gravitation corrections are not small 
enough to suppress it adequately. 

If 1/Mp terms are forbidden (for example by some unbroken symmetry), then 



1/Mp terms (d=6 operators) become important. An example of such operator is 



This term gives an e which is about 10 orders of magnitude less than needed to 
produce observable effects. 

Therefore, while in the MSSM l/M Pl terms are too small, the 1/Mp; terms are 
too large and need additional suppression, i.e. small (3. 

If wormhole effects are responsible for the terms we are discussing, they can 
contain a topological suppression leading to very small (3. Generically, this sup- 
pression is described by an e~ s factor, where S is an action of a wormhole which 
absorbs the Rp charge. In the semi-classical approach S ~ 10. In the case of the 
Peccei-Quinn symmetry such estimates give S ~ \n(M P / f PQ ) « 16 resulting in a 
suppression factor j3 ~ 10 7 , which is needed in our case. A detailed discussion 
of such suppression factors in wormhole effects is given in [|T7[]. It is shown that 
the action S is connected with the size of the wormhole throat R(0) and can vary 
from S « 6.7 for the naive estimate R(0) « Mp 1 , up to Six 2 /g 2 str « 190 in string in- 
spired models. Thus if the action is close to its semi-classical value the topological 
wormhole effects can provide the suppression needed for the long-lived neutralino 
13 ~ 10~ 7 — 10~ 5 to satisfy the astrophysical constraints. 

Apart from topological wormhole effects, suppression of d=5 operators can occur 
due to some additional symmetry. Let us assume that there exists a singlet sec- 
tor which communicates with the MSSM sector only gravitationally through non- 
renormalizable terms in the Lagrangian. R-parity can be broken spontaneously in 
this sector, for example, due to some i?p-odd field r] developing a non-zero VEV. Rp 
violation can penetrate the MSSM sector through non-renormalizable interactions 
between r] and the MSSM fields, n contrast to the first case, gravity is not directly 
responsible for the breaking of Rp but it leads to effective Rp violation in the ob- 
servable sector through the presence of non-renormalizable interactions. We now 
discuss this possibility of hidden R-parity violation in a model independent way 
and then provide an example. 

Let us assume the existence of an SU(3) x 577(2) x [7(1) singlet field 77 coupling 
to the MSSM fields only through non-renormalizable terms. This can be achieved 
by a proper symmetry as we shall discuss. There are four dimension-five operators 
involving rj which lead to R P violation: 



M P 



{{LH 2 ){H x H 2 f) 



D 



(10) 



O x = §;{U c D*D c ri) F 
3 = $-(QD°H lV ) F 



2 = ^- p {LLE% F 
4 = ^-(LH 2V *) D 



(11) 



where 01,2,3,4 are parameters of order one. These operators conserve R P if the field 
77 is chosen odd. The non-zero VEV (77) then breaks R P and leads to the effective 
interactions in eq.(2) with couplings given by (i =1,2,3): 

\ = a i (v)/Mp e = a 4 (F v *)/M P (12) 

The effective Rp violation among the MSSM fields is governed by two distinct 
scales. One scale is (77) and determines the trilinear interactions of eq.(2), while the 
other scale (F v * ) characterizes SUSY breaking in the singlet sector and determines 
the bilinear term e. In general, these two scales could be quite different. The 
constraints derived in ((|8|)) imply 

(77) < HrtJeV; (Fr) < 10 2 GeV 2 (13) 

If SUSY is broken through the usual soft terms, the constraint on (F v ) can easily 
be satisfied, as we will show in a specific example. In contrast, the constraint on 
the trilinear coupling implies very small 77 VEV which may be unnatural. To avoid 
this one should forbid the corresponding dimension-5 operators, in which case the 
dominant Rp violation would arise from dimension-six interactions. For example, 
the operator 

-^[{LH 2 ){H x H 2 )rj\ F (14) 

leads to e ~ 10 34 (77) GeV The effect of this term could be observable only if R P 
violation in the singlet sector occurs close to the unification scale. 

Let us now consider a concrete realization of the above scenario. We will impose 
a gauged discrete symmetry in order to prevent gravitational breaking. -Rp by 
itself can be gauged in a suitable extension of the MSSM but it cannot prevent the 
occurrence of a term like LH%r\ which lead to a large Rp violation when 77 develops 
a vev. We therefore look for a more general symmetry which leads to effective 
i?P conservation at the renormalizable level. Gauged Z N symmetry of the type 
considered in [^] provides an example. The 6 is assumed to carry the Z^-charge 
-1. The Z N charge of one of the observable superfields is set to zero by appropriate 
redefinition of the Z N generators. Requiring that the standard R P conserving 
couplings of the MSSM fields are allowed by the Z N symmetry we determine the 
charges of the remaining fields in terms of two parameters x and y according to: 

Q U C ,H! D C ,H 2 L W_ I 3 
x 2-xy2-(x + y)2f 



where we have introduced a singlet field Y in addition to rj, in order to obtain R P 
violation. With these charge assignments dimension-4 terms respect R P and the 
r], Y do not couple to the MSSM fields in the renormalizable Lagrangian, as long 
as x 7^ 2 and x — y ^ 0, —2, N/2. The most general Z^v-invariant renormalizable 
superpotential in this case can be written as QI5|] 



W = W MSS M + SY(r] 2 -f 2 ) (15) 

The above superpotential gives rise to a VEV for r] at the supersymmetric mini- 
mum, leading to effective R P violation for the MSSM fields through operators of 
dimensionality > 5. The choice 2 + y — x = N/2 allows the dimension-6 operator of 
(|T3]). The allowed higher dimensional terms are given in this case by 

Cnr = -j^(LH 2V *) D + ±[5i(LLE°rr)D + 
5 2 (QD c Lr]*) D + 5 3 (LH 2V *Y) D + 5±{LH 2 H x H 2 ri) F ) (16) 

Note that the dimension-5 operator above cannot be forbidden if the dimension-6 
term in (|TJ) is to be allowed. As discussed above it does not lead to large R P violation 
as long as SUSY remains unbroken in the singlet sector, as happens with the 



superpotential choice in fll5|). In a realistic situation soft SUSY breaking introduces 
terms which make F Vt y non-zero. If (rj) < M S usy ~ TeV then (F v ) ~ 2 ^/ 2j4 where 
my ~ A characterizes soft SUSY breaking. Choosing 5 ~ 10~ 2 , my ~ 10 3 GeV, 
A ~ 10 2 GeV and (rf) < 100 GeV we have F v < 10 2 GeV 2 in agreement with the 
constraint of (|T3|). If (rj) is much larger than M S usy then (F v ) would also be large 
and induce large e. This can be prevented by a symmetry. Specifically, if the kinetic 
energy terms for the singlet fields are chosen to be no-scale type [ ]2~0"[ ] then (F^y) 
vanish at the minimum of the potential and effective breaking of R P would arise 
only from the dimension-6 operator. This operator could lead to observed signatures 
if (77) is very large, near the unification scale. 

The Z N symmetry introduced above can be gauged if it satisfies discrete gauge 
anomaly constraints discussed in p(J[j. In our case these are given as: 



-2N g + 6 = h N (17) 
N g (y-4)+4 = k 2 N 
N g (-7 + y-x) + N/2 - 9 = k 3 N + k k 4 N/2 



where k is 1(0) for even (odd) N and £4,2,3,4 are integers. The first constraint is 
automatically satisfied for three (N g = 3) generations. The remaining constraints 



can also be satisfied for appropriate choices of x, y and N. A choice which satisfies 
all the anomaly constraints and leads to the required interactions in Eq. ( |I5| ) and 
is given byiV = 3,x = l/6 and y = —1/3. Clearly many other choices would be 
possible. 

In summary there is no deep theoretical motivation for R-parity to be abso- 
lutely conserved. In this paper we have analyzed the possibility that R-parity can 
be weakly violated by gravitationally induced terms suppressed by inverse powers 
of the Planck mass. Due to astrophysical constraints (mainly to the positron pro- 
duction in our Galaxy) this extremely weak R-parity violation is compatible with 
hypothesis that the neutralino is a DM particle only if its lifetime is about 10 10 
times longer than the age of the universe. We analyzed gravitationally induced 
dimension-5 operators and demonstrated that they break R-parity too strongly to 
comply with this constraint. We showed that the unstable neutralino as a DM 
particle is possible if dangerous d=5 terms are strongly suppressed. For the MSSM 
case, the required suppression could be provided by the classical wormhole action 
7. 

Another possibility for very weak R-parity breaking can be provided by existence 
of additional gauge discrete symmetry. We constructed a model with a gauge Z N 
symmetry and SU(3) x SU(2) x U(l) singlet field rj, which communicates with the 
MSSM fields only through gravity. R-parity violation is driven by VEV of this field. 

The decaying neutralino has interesting astrophysical signatures. In some mod- 
els [ {Tip the neutralino decay to the Majoron J, \ —> ^+ J may be dominant, resulting 
in a detectable isotropic flux of mono-energetic neutrinos. In the more general case 
of R P breaking by d=5 operators discussed above, the neutralino decay signature is 
weaker and is given by the ratio of the signals from the Sun and Earth to that from 
the Galactic halo. The signal from annihilation of neutralinos in the Earth and the 
Sun is the same as for a stable neutralino, while the positron and antiproton fluxes 
from the Galactic halo could be strongly enhanced due to neutralino decay. 
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